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SUMMARY 

Magnetic field intensity and currents passing through the coils of the National Aerospace 
Laboratory (NAL) 10cm Magnetic Suspension and Balance System (MSBS) were measured while 
a cylindrical model was oscillated along x,y,z and also about y and z axes, respectively. The 
model was made of alnico 5 and was 8mm in diameter and 60mm long. Two kinds of tests were 
carried out. Amplitude of the oscillation was varied at a frequency of 10Hz. Frequency was 
varied from 1 to 50Hz in the other test. Results of the tests show that the relation between coil 
currents and magnetic force acting on the model is affected by frequency. They also show that the 
relation between measured magnetic field intensity and the force in vertical direction is 
independent of the frequency below 30Hz. Using the measured magnetic field intensity, the 
vertical force can be evaluated at the MSBS instantaneously when a model moves at frequencies 
below 30Hz. A static drag force calibration test was carried out at the 60cm MSBS. Obtained 
relationships between measured drag coil currents and loads shows large hysteresis. 

INTRODUCTION 

Magnetic Suspension and Balance Systems (MSBS) are a kind of model support system of 
wind tunnels for supporting a model in flow with magnetic force. It can avoid the model support 
interference because the flow field in a test section will not be affected by the magnetic field except 
for some very high speed flow, etc. Magnetic field for suspending a model is generated by 
currents passing through some coils arranged around a test section. Then the magnetic force acting 
on the model corresponds to the currents. When the model is at a fixed position in the flow, the 
magnetic force must balance with the aerodynamic force acting on the model plus the gravity 
force. This means that the system can work as a balance for the aerodynamic force by measuring 
the coil currents. The relation between the magnetic force and coil currents can be decided 
uniquely by some calibration tests (static force calibration tests) in almost all cases. In case of a 
model in motion, the difference between the magnetic force and the aerodynamic force places the 
model in motion. Then the aerodynamic force can be evaluated by subtracting the force driving the 
model motion from the magnetic force. The inertia force can be estimated by measuring model 
motion. This means that the system can measure dynamic force. It is also easy for the system to 
get accurate model position because it uses the position data for its control. It is easy to create 
suitable forces on a model by placing the model in suitable motion in the system. The tests 
corresponding to the magnetic force to coil currents are called dynamic force calibration tests. 
Dynamic calibration tests can replace static force tests because there is no difference between the 
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force driving the model motion and the force by pulling with threads in the sense of the force 
acting on the model. Some dynamic force calibration tests were carried out at the 10cm MSBS. 
Preliminary static force calibration tests were conducted at the 60cm MSBS. 

SYMBOLS 

g gravity acceleration ( = 9.8m/sec 2 ) 

H magnetic field intensity, (H x ,H y ,H z ) (T) 

Hi a magnetic field intensity component at a position 
along the coil no.i axis. (T) See Figure 2. 

Hmonitor x component of the magnetic field intensity inside test section (T) 

Hdrag, Hside, Hlift, Hpitch, Hyaw See equation (1) 

Ii current passing through coil no.j (A) 

Idrag, Iside, Ilift, Ipitch, Iyaw See equation (1) 

Iyy, Lz moment of inertia about the y and z axis, respectively, (kgm 2 ) 

Kdrag, Kside, Kiift, Kpitch, Kyaw proportional constant between force and magnetic 

field intensity combination. See equations (2) and (3). 

m model mass (kg) 

nix x component of magnetic moment (Wbm) 

(x,y,z) coordinate system. See Figure 1. 

6 angle about the y axis. 

<]) angle about the z axis 

DYNAMIC FORCE CALIBRATION TEST 
Experimental Design 

Magnetic field around the test section together with coil currents was measured with Hall 
sensors during dynamic force calibration tests. Some preliminary test results have been published. 
As mentioned in References 1 and 2, even the most suitable combination of measured coil currents 
showed differences from model position in the sense of phase. The difference was observed to 
depend on the motion frequency. On the contrary, a suitable combination of measured magnetic 
field intensity components appeared completely in phase with the model position. In this 
dynamic calibration test, forced oscillations of a model were performed in two ways. Frequency 
varied from 1 to 50Hz along the x, y and z axes and about the y and z axes in one way. The 
amplitude of 10Hz oscillation varies in the other way. In order to avoid the effect of eddy 
currents on metal test section walls, the test section was removed during the test. 

The model is 60mm long and 8mm in diameter and of alnico 5 permanent magnet. In order to 
measure the model position, the model was wrapped with thin white paper. The paper measured 
0.2g +/- O.lg in mass. The shape is 60mm high x 52mm wide x 0.1mm thick. A 4mm wide black 
line was printed along the center line of the paper. The line is used to measure the x position of 
the model. The whole mass of the model including the paper was measured as 23. 8g +/- 0.1 g with 
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a balance. The inertia moment of the model was evaluated at 2.87 x 10 5 kgm 2 about an axis through 
the gravity center normal to the model axis by calculation from the model shape and mass. The 
model was supported by the control of 5 degrees of freedom except for rolling motion about the 
model axis. Besides, constant currents of 4 A in magnitude pass through the 4 side coils to 
generate an additional constant magnetic field which makes model position stable in y and <{>. 

The model was suspended at the center of the 10cm MSBS which is the coordinate system 
origin. All oscillations were pure sinusoidal, pure heaving and pure pitching motions. In the 
pitching (and yawing) motions, the center of gravity was controlled to be kept at the origin. 

Coil currents were measured from monitor outputs of 10 power amplifier units. In order to 
evaluate the monitor output accuracy, they were compared between 0 to 50Hz with the currents 
measured with a zero inductance Shunt type resistance of 0.1 Q and of 0.1% accuracy. The 
monitor outputs are less than the coil currents by about 1.2% but the difference is independent of 
the frequency. The monitor outputs were delayed against the real coil currents by 0.5 degree at 
the highest frequency of 50Hz. 

The monitor outputs were measured and recorded in a personal computer with 12 bit AD 
converters with sample and hold function which are 0.2% accurate over a +/- 10V range. All 
converters were adjusted to 5mV error at most in their whole range with a reference voltage 
generator of 0.01% accuracy. The coil currents and the magnetic field intensity components were 
measured at the same time. 

The Hall sensors, THC126 (Toshiba), were used to measure the magnetic field intensity 
components at the 8 points shown in Figure 2. H x is symbolized as H 9 here. Although the 
measurement of H 9 is not disturbance-free against the flow field, it is possible to measure H 9 
like a wake survey in wind tunnel tests without serious affect on the test result. The sensors were 
driven in the constant current mode. They were calibrated with the model 9903 of F. W. BELL 
which is 0.1%FSR accurate. They showed good linearity between the Hall sensor outputs and the 
Gauss meter output. An example of the calibration test results are shown in Figure 3. 

The model position was measured with an in-house-developed model position sensing system. 
The system was described in detail in Reference 3. The position sensor was calibrated by a 
calibration model of 8mm diameter wrapped with the same paper as the dynamic force calibration 
model. It was positioned on a stage which can very all positions in 6 degree of freedom. Some 
examples of the position sensor calibration test are shown in Figure 4. There are some 
interferences between measured positions, particularly in y and <|>. The interference will induce 
some unexpected motion of the model. 

The control speed of the MSBS measured 489.8Hz. The obtained model position by the 
sensor is delayed by 3.0ms from the measuring time of the coil currents because it is estimated 
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with CCD sensors as presented in Figure 5. The times and frequencies shown in the Figure were 
measured with a universal pulse counter. 


Dynamic Force Calibration Analysis 


The model is assumed a perfect cylinder and is also assumed to be magnetized along the model 
axis like (m x , 0, 0). The magnetic force components along the x, y and z directions acting on the 
suspended model can be approximated by the magnetic field intensity gradient about the model 
multiplied with m x The moments about the y and z axes can be also approximated by the 
magnetic field intensity about the model multiplied with (m x , 0, 0). Then the moment components 
about the y and z axes are approximately proportional to the averaged H z and H y over the 
model, respectively. To monitor the magnetic field intensity about the model, the following 
combinations of the measured magnetic field intensity components are defined: 


Hd, 


H 1 +H 2 -H 3 -H 4 +H 5 +H 6 -H 1 - Hz 


'ratf= 


8 


Hmonilor — 7 j 9 , 


Idrag — 


7o +h 


IT H 2 + H 4 — H 6 — H 8 , Il + H+H+Is 

/j side = , l side — 
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Hi, f t = 


Hi + Hi- -Hi- Hi 


I hft = 


h + h + Is + I- 


H pitch ~ 


Hx + Hi + Hi + H- 


Ipitch — 


I\ + 1 31 — 1$ — 1 7 


Hvcrw — 


Hi + H 4 + H 6 + Hs 


lyaw — 


1 2 + 1 4 — If, — Is 


The motion of the model gravity center satisfies the following equations: 

mk = Kdrag' Hdrag, Of" fflX — Kmonitor' H monitor, 

rriy = Kstde • Hs,de, ... ( 2 ) 

mz = Kiifi • Hitft - mg, 


where K’s are proportional constants between the quantities defined above and the corresponding 
force. The motion of model rotation about the y and z axes at the gravity center satisfies the 
following equations: 


Iyy 0 = K pitch ' H pilch , 
Izz iji = IC yaw ’ H yaw , 


(3) 


where K’s are the similar ones as above mentioned. If motion is sinusoidal, the force driving the 
model motion varies sinusoidally as the motion. Besides, the motion and the force must be in 
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phase. The model was forced to make as pure a sinusoidal motion as possible in the dynamic 
force calibration test. By examining the phase between the motion and the quantities defined in 
the above equations, it can be estimated whether or not the quantities are suitable for evaluating 
the dynamic force. 


Dynamic Force Calibration Test Results 


Heaving Motion 

Figure 6 shows trajectories of H^, Inn and z position of the model with respect to time 
during a 10Hz heaving motion. The three quantities are normalized with their rms values. H llft 
and z are perfectly in phase. The waveforms of the two also look like a single frequency 
sinusoidal oscillation. Figure 7 shows also H ljft , Ii ift and z vs. time during a 20Hz heaving 
motion. The same observation can be remarked as in the 10Hz motion. This fact confirms that 
Hjift is directly proportional to the force in the z direction. On the contrary, there is observed 
apparent phase difference between -Inf t and z. This suggests that 1^ is not proportional to the 
force unlike Hi ift . 

The rms value of the force during the heaving motion was evaluated from the measured z 
position of the model. Figure 8 shows the relation between Hiift and the force in their rms values. 
The symbol of x shows the case of changing amplitude of 10Hz heaving motion. A dotted line in 
the figure is a least square approximation line fitted to the results in the amplitude change case. 

The symbol of an open circle shows the case of changing frequency between 1 to 30Hz with 
various amplitudes. The open circles are on the approximation line except for frequencies less 
than 5Hz. The results at frequencies higher than 30Hz are not on the line. One of the causes is 
poor magnetic field control because of large induced electromotive force. The maximum error of 
the dynamic force balance in the z direction reduces to about 2.4% in the region of the figure. 

The error is observed in a 1Hz heaving motion and its value is 2mN. The change of coil currents is 
about the resolution of the current amplifier units at the frequency. Then, suitable control of the 
magnetic field cannot be expected. The fact suggests that a core permanent magnet must be chosen 
to meet test requirements. When a small force in the z direction is measured, a weaker and light 
magnet must be used. 

Figure 9 also shows that relation between the force and I Uft in their rms values. The symbols 
are the same as in Figure 8. Although the results in the case of amplitude change are on a line, the 
results in the case of frequency change are not. They show a frequency dependence when the 
frequency is higher than 10Hz. It means that the obtained data must be compensated by the effect 
of frequency dependence if dynamic force is estimated with I^ft. The magnetic circuits of the 
10cm MSBS are made of iron blocks and have hysteresis loss dependent on frequency. This is 
suspected to be the cause of the frequency dependence. 
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Pitching Motion 


Figure 10 shows trajectories of H pilch , l pjlch and 0 of the model attitude with respect to time 
during 10Hz pitching motion. The three quantities are normalized with their rms values. H P i tdl 
and 0 are in phase. The waveforms of the two also look like single frequency sinusoidal 
oscillations. Figure 1 1 shows also H pitch , I pitch and 0 vs. time during 20Hz pitching motion. 
Although H pitch and 0 are in phase independently of frequency, H pitch trajectory differs from a 
pure sinusoidal waveform around its peaks. H pitch is approximately proportional to the torque 
about the y axis according to the figures. On the contrary, there is observed apparent phase 
difference between I pitch and 0. This suggests that I pitch is not proportional to the torque unlike 

Hpitch- 

The rms value of torque about the y axis (pitching moment) was evaluated from the measured 
0. Figure 12 shows the relation between the pitching moment and H pitch in their rms values. The 
symbols are the same as in the case of the heaving motion test. The open circles are not on the 
approximation line. The similar difference between the two kinds of tests is found in the relation 
between pitching moment and Ip itc h as shown in Figure 13. The difference is larger than in 
Figure 12 with pitching moments larger than 0.7mNm. It suggests that the frequency dependence 
lies in the moment evaluation by l pitc h as in the heaving motion test. The open circles and cross 
symbols in the range less than the pitching moment of 0.2mNm are not on a line. The cause is that 
Ipitd, and H pi , ch are very small and nearly out of the controllable range. Figure 14 shows the 
Hpitch, Ipitch and 0 vs. time during 4Hz pitching motion. The model position is controlled but 
Hpj te h and Ipitch are scattered around the pitch angle trajectory. The other open circles are 
approximately on a line but the line is different from the line approximated by the cross symbols. 
It means the relation between pitch angle and H pitch depends on the frequency. The cause of it 
has not yet been isolated. 

Oscillation in the x direction 

Figure 15 shows H^g, I drag, and H monitor and x position with respect to time during 10Hz 
oscillating motion in the x direction. The four quantities are normalized with their rms values. 
Either of four is not perfectly in phase with z. H mon ; t or is the most in phase with z of the three. 
The rms value of the force in the x direction was evaluated from the measured x position of the 
model. Figures 16, 17, and 18 show the relations between the force and H^g, Idrag, H m omtor> in 
their rms values. Results in the frequency change are nearly on a line approximated by those in 
amplitude change at 10Hz in the case of using H mo „i, 0 r- The other two cases show the frequency 
dependency. But results in the case of frequency change are nearly on a line approximated by 
those in amplitude change only around 10Hz in case of using Idrag- One of the causes is the air- 
cored drag coils without the hysteresis loss by iron cores. 
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Oscillations in the y direction and about the z axis 

Figure 19 shows H side , I side and y position of the model with respect to time during 10Hz 
oscillation in the y direction. The three quantities are normalized with their rms values. Three 
waveforms do not look sinusoidal at all. It means that the control is not adequate. Cause of the 
distorted waveforms is the poor accuracy in y of the model position sensor. The rms value of 
force in y was evaluated from the measured y position of the model. Figures 20 and 21 show 
the relations between the force in the y direction and H sidei I s ; de in their rms values. Results in 
the two test cases are scattered around a line which is approximated by the amplitude change tests 
at 10Hz. 

Figure 22 shows H yaw , I yaw and <|> position of the model with respect to time during 10Hz 
oscillating motion about the z axis. The three quantities are normalized with their rms values. 
Three waveforms do not look sinusoidal. It means that the control is not adequate. Cause of the 
distorted waveforms is the poor accuracy in <}) of the model position sensor. The rms value of 
yawing moment was evaluated from the measured (|> of the model. Figures 23 and 24 show the 
relations between the yawing moment and H yaw , I yaw in their rms values. Results in the two test 
cases are scattered around a line which is approximated by amplitude change tests at 10Hz. 

STATIC FORCE CALIBRATION TEST AT THE 60CM MSBS 

A preliminary static force calibration test was carried out at the 60cm MSBS. The 60cm 
MSBS is the largest in its test section size, which is described in References 1 and 2. It is 
operated in 5 degree of freedom control. But y position and yaw angle control is very slow 
because of poor accuracy of the model position sensor in those directions. Only axial force was 
calibrated by pulling the model with lOOg scale weights and measuring the drag coil currents. 
Figure 25 shows a picture of the test. The model is 381mm long and 55mm diameter and it 
contains a cylindrical permanent magnet core of Fe-Cr-Co magnet which is 50mm diameter and 
300mm long. In order to measure the model position, the model was wrapped with a thin white 
paper as in the case of the 10cm MSBS model. The whole mass of the model measured 5.1kg. 

The model was supported by the control of 5 degree of freedom except for rolling motion. Drag 
coil current was measured with monitor output of a power amplifier unit for the drag coil. The 
monitor outputs of the 5 power amplifiers are accurate of 0.1% FSR up to 75A for constant 
output. The monitor output was measured and recorded by a FFT analyzer. The model position 
was measured with the same typed model position sensing systems at the 10cm MSBS. The 
position sensor was calibrated with the model. It was positioned on a stage which can very all 
positions in 6 degree of freedom. The control speed of the MSBS measured 248Hz. The obtained 
relation between the drag coil current and applied load is shown in Figure 26. Apparent hysteresis 
is observed. Drag force calibration tests were carried out several times. Hysteresis is observed 
every time. Figure 27 shows the results of another test in which the axial load increases and 
decreases in small steps of about lOg. Large hysteresis is observed. 
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REMARKS 


Magnetic field intensity and coil currents of the 10cm MSBS were measured when a 
cylindrical model was oscillated along x, y, z and also about y and z axes, respectively. Two 
kinds of tests were carried out. Amplitude of the oscillation varies with constant frequency of 
10Hz in one mode. Frequency varies from 1 to 50Hz in the other mode. 

In case of heaving motion, H uft and z position are in phase but Ii ift and z are out of phase. 
The relation between H !in and the force in z direction is very linear in the tested range and is also 
independent of the frequency below 30Hz at the 10cm MSBS. On the contrary, the relation 
between I| ift and the force is not linear but also depends on the oscillation frequency. H| jft can be 
used for estimating lift force acting on a model in motion. 

In case of pitching motion, 0 trajectory looks like a pure sinusoidal waveform but H p j tc h 
trajectory differs from a pure sinusoidal wavefonn around its peaks although H p j tc h and 0 are in 
phase. I p i tc h and 0 are out of phase. The relation between H pitch and pitching moment is linear 
but depends on the test modes. The cause of the test mode dependence has not been isolated yet. 
The relation between Ii ift and the moment is not linear but also depends on the frequency. 

x position trajectory looks like a good sinusoidal wavefonn but H^g, H moi iitor and Idrag 
trajectories do not look like good sinusoidal wavefonns. x and either of the three is not 
apparently in phase. The relation between H^g and the force in the x direction is linear but 
depends on the two test modes. H mon j tor and the force is linear and is also independent of the test 
mode. It suggests that the relation is independent of the frequency. The relations between I drag 
and the force and between H^g and the force are linear but depend on the test modes. 

Oscillations in the y direction and about the z axis are not pure sinusoidal ones because of 
poor accuracy in y and <J> of the model position sensing device at the 10cm MSBS. The 
relations between H sjde and the force in the y direction and between I side and the force look 
linear but are not reliable to estimate the force from the measured H sid e or I si d e . The relations 
between H yaw and the yawing moment and between I yaw and the moment also look linear but are 
not reliable either. 

At the 60cm MSBS, drag force calibration tests were carried out. The results show large 
hysteresis between load and drag coil current. The cause has not been isolated yet. 
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Hall sensor 

Figure 2 Hall sensor position 


Hall element No. 1 y = 0. 04x - 0. 0062 



Gaussmeter output (G) 


Figure 3 Hall sensor calibration test result 
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Figure 4 model position sensing device calibration test result 
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Figure 5 data acquisition timing diagram 
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Figure 25 drag force calibration test 
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